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Summary 


The purpose of this thesis is to analyze the divergent behaviour of perturbative series in QCD, 
in order to estimate the uncertainty on finite-order perturbative results. With a new LHC 
era starting this year (2015), precision in theoretical predictions is important more than ever, 
especially for the newly discovered Higgs boson. Unluckily, the perturbative series of the cross 
section of Higgs production in gluon fusion (the most relevant channel at the LHC) converges 
rather slowly. Furthermore, the methods we have to determine the size of higher order corrections 
have a poor accuracy, and most of them lack of a physical interpretation. 

Our starting point is that any observable in QCD or QED which is a perturbative expansion 
in the coupling constant a* or a is divergent, and this has been known since a long time (Dyson, 
1952). In Chapter we will first illustrate Dyson’s argument. Then, we will introduce some 
tools that can help us deal with divergent series, in particular Borel summation, and see some 
instructive examples. Finally, we will talk about the three most important known sources of 
divergence in QCD: the Landau pole, renormalons and instantons. 

The Landau pole divergence is closely related to the concept of resummation, which we shall 
introduce in Chapter A certain class of corrections to the leading order cross section, the 
emission of soft gluons, is known to produce logarithmically enhanced contributions, which can 
be resummed to all orders. However, this sum is convergent in Mellin space, but cannot be 
transformed back in physical space, due to a branch cut for large values of the Mellin moment 
N. The presence of the Landau pole is ultimately responsible for this branch cut. As a result, 
the expansion of the resummed cross section in physical space, which happens to be (at least for 
the Higgs case) a very good approximation to the exact hxed-order calculation, is a divergent 
series. 

In Chapter we will see that another class of perturbative corrections, the fermion bubble 
diagrams, produce a factorial divergence when integrating over very high (UV) or very low (IR) 
internal momentum. This divergence is called renormalon divergence and we will see how one 
can generalize the concept of renormalon, dehning it as a pole in the real axis of the Borel plane, 
located at an integer multiple of /3o. 

There are other poles in the Borel integration path: the instantons. They produce another 
factorial divergence which is this time related to the growing number of Feynman diagrams 
with increasing perturbative order. In some quantum mechanical systems and even in some 
quantum field theories, one can sum the instanton-anti instanton contribution to the divergent 
perturbative expansion and obtain a so-called resurgent trans-series, eliminating completely the 
divergence. Unluckily, this is not the case of QCD, where IR renormalons are much closer to 
zero (and therefore lead to a faster divergence) than instantons. 

After introducing the sources of divergence, in Chapter we will illustrate the current most 
used models for theoretical uncertainties. Since most of the theoretically calculable cross sections 
are known up to the third or fourth perturbative order at best, those methods make an attempt 
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to describe the effect of the unknown higher orders. An issue which we will not discuss in 
detail but deserves a special mention is the impact of parton distribution functions (PDFs) on 
theoretical uncertainties. 

Finally, in Chapter we will derive our own model for estimating theoretical uncertainties, 
based on what we know about the sources of divergence mentioned above. We will then apply 
that model to two of the most relevant processes at the LHC: Higgs production and tt produc¬ 
tion. Our first concern will be to compute the perturbative order at which each of the sources 
of divergence starts to kick in. After that, we will compute our estimate for the theoretical 
uncertainties on the last known order and we will compare our result with that of the already 
known models. 
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Chapter 1 

Divergent series in QCD 


The Large Hadron Collider has just restarted (early 2015) after two years of maintenance and 
upgrading. As we write, the center of mass energy of the collisions is set to 13 TeV, which means 
that the cross sections of most of the events that we can observe are significantly increasing with 
respect to a couple of years ago. With no evidence of Beyond the Standard Model physics at 
8 TeV, precision measurements in the SM are needed from an experimental and a theoretical 
point of view. 

Let us focus on theoretical accuracy, starting with an example. The cross section for the 
production of a Higgs boson in gluon fusion at LHC @13 TeV is, as any observable computed 
in quantum field theory, a perturbative series. This means that if we want to give a theoretical 
prediction for the cross section, its form will be 

agg^H = Co + cias + C2al + ..., ( 1 . 0 . 1 ) 


where Ug is the strong coupling constant, which computed at the Higgs mass scale has the 
value asirnn) = 0.1126. Being ag a perturbative parameter, one could expect that the difference 
between, for example, the Next to Leading Order (NLO) and the NNLO should be small. If this 
were case, the calculation of high perturbative orders, which usually requires a lot of time and 
effort, would be only needed if one wanted to have an extraordinary precision. However, as we 
can see in Table [lH this perturbative series is slowly convergent at the first known orders. The 
perturbative corrections can even be bigger than the whole cross section at the previous order. 
The situation does not improve changing the energy of the collision. This is one issue we want 
to address: if we want to give a good theoretical prediction to the Higgs cross section, we need 
to take into account in some ways even the unknown higher orders. 

Another important issue is the fact that series like (1.0.1) are known to be slowly convergent. 
Actually, Dyson in 1952 showed that any perturbative expansion in QCD is divergent. As a 
matter of fact, we know well some physical sources of divergence: instantons, renormalons and 
the Landau pole in soft gluon resummation. What we do not know is the perturbative order 
at which these divergences start to kick in. If, for example, the series (1.0.1) started to diverge 


c.m.e LO [pb] NLO [pb] NNLO [pb] % 

8 TeV 5.37 12.78 17.15 138.18 % 34.21 % 

13 TeV 12.25 29.40 39.24 139.98 % 33.46 % 


Table 1.1: Higgs cross section in gluon fusion at 8 TeV and at 13 TeV 
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at the 4**^ perturbative order, then any exact N^LO calculation would not make sense, because 
that result would be distant from the exact value of the series. 

This chapter is so organised: first we will present Dyson’s argument and show why all 
perturbative series in QED and QCD are divergent. Then we will describe the basic tools we 
can use to treat divergent series. 

1.1 Dyson’s argument 

In 1952, Freeman Dyson presented a very straightforward argument which led to the conclu¬ 
sion that all the power series expansions in use in quantum electrodynamics are divergent after 
the renormalization of mass and charge. This same argument can be applied with just a little 
effort to QCD. 

The starting point is that all existing methods of handling problems in quantum electro¬ 
dynamics give results in the form of power series in a = e^/(47r). Let us consider a generic 
observable a: 

a{e^) = ttQ + a2e^ + + - (l-l-l) 

Of course, the coefficients Oj are finite after mass and charge renormalization. The series as a 
whole, though, cannot be treated with the same techniques. If the series converges, its sum is a 
calculable physical quantity. But if the series diverges, it becomes difficult to calculate or even 
of define the quantity which is supposed to be represented by the series. 

We know that QED is equivalent to a theory of the motion of charges acting on each other 
by a direct action at a distance. The interaction between two like charges is proportional to e^. 
Suppose now that the series (T(e^) converges for some positive value of e^, this implies that (T(e^) 
is an analytic function of e at e = 0. Then, for small values of e, cr(—e^) will also be an analytic 
function with a convergent power series expansion. 

But we can also find a physical interpretation for cj(— e^): it is the value that would be 
obtained for a if the interaction between charges of the same type had a minus sign. In a ficti¬ 
tious world like that, charges of the same type attract each other and the classical macroscopic 
potential is just the Coulomb potential with the sign reversed. But in these conditions the vac¬ 
uum state is not the state of lowest energy. In fact, one could construct a pathological state by 
creating a large number of electron-positron pairs and bringing the electrons and the positrons 
in two separate regions. In a state thus made, the negative potential energy of the Coulomb 
forces is greater than the sum of the total rest energy and kinetic energy of the particles. 

This can be done without using particularly small regions or high charge densities, so that 
the validity of the classical Coulomb potential is not in doubt. Let us suppose now that a system 
is given at a certain time with only a few particles present. There exists a high potential barrier 
separating this physical state from the pathological state described above. However, due to 
the quantum mechanical tunnel effect, there is a finite probability that the system will evolve 
towards the pathological state. Therefore, any physical state is unstable against the spontaneous 
creation of infinite particles. Furthermore, once a system finds itself in a pathological state, there 
will be an inevitable creation of more and more particles. In these conditions it is impossible 
that the integration of the equations of motion, starting from a given state of the fictitious world, 
should lead to well defined analytic functions. Therefore a{—e^) cannot be analytic and the the 
series cannot be convergent. 
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1.2 Divergent series 


We have just seen that any perturbative series in QED or QCD is intrinsically divergent. Even 
so, there is much that we can say about divergent series and the interpretation of their sum. 
An exhaustive review about the treatment of divergent series can be found in Ref. [^, we will 
mainly follow Appendix D of Ref. for our purposes. Let’s start with some definitions. A 
generic series 

OO 

S = J2ck ( 1 . 2 . 1 ) 

k=0 

is convergent if, being its partial sums 

n 

Sn = ^Cfc, (1.2.2) 

k 


the following limit 

s = lim Sn (1.2.3) 

n^oo 

is finite. In this case such limit s is called the sum of the series and we say that S = s. Otherwise, 
we say that the series is divergent. Another very useful dehnition is that of absolute convergence: 
a series is absolutely convergent if the series of the absolute values \ ck\ is convergent. 

If the series is a power series, 

S{z) = '^Ckz’^, (1.2.4) 

k 

then its convergence depends on the value of In particular, what happens is that the series is 
convergent for values of z that lie in a circle of radius r, where 


r = lim 
k—^oo 


Ck 

Cfc+1 


(1.2.5) 


The convergence of the series implies the analyticity of the sum s(z) inside such circle. Therefore, 
if we expand a function f(z) around some point zq, the radius of convergence of the expansion 
can be at most the distance between zq and the singularity that is closest to it. 


1.2.1 Asymptotic expansions 

Sometimes, a function f{z) might admit a series expansion around a certain value of z that is a 
divergent series. In other words, if we expand f{z) around z = 0, we have that a series expansion 


Siz) = '^CkZ^ 
k 

is asymptotic to f{z) if there exists a constant K such that 


\f{z) - s„( 2 :)| < K Cfe+i 


( 1 . 2 . 6 ) 


(1.2.7) 


for all n. 
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1.2.2 Sum of divergent series 

Whenever we have a divergent series, by definition the limit of the partial sum is infinite (or 
does not exists). However, this has nothing to do with the finiteness of the sum of the series. 
In fact, we can define the sum of a series in another way, it doesn’t have to be the limit of the 
sequence of the partial sums. 

As an example, let’s consider the following (divergent) series 


P = ^(-1)'’’= 1 - 1 + 1 - 1 + .... 


k=0 


This series is divergent because the partial sums 


Pn = 




k=0 


0 for even n 
1 for odd n 


( 1 . 2 . 8 ) 


(1.2.9) 


oscillate between 0 and 1, and therefore the limit for n —)• oo ofpn is not defined. However, there 
are multiple arguments that assign to the sum the value 1/2. For example, we can manipulate 
the definition of the series to obtain an equation for P: 


P = ^(-1)^ = 1 + J^(-l)^ = 1 + 1 = l-P (1.2.10) 

fc=0 k=l k=0 k=0 

from which we get P = 1/2. Another way to obtain P = 1/2 is to consider the power series 


Piz) = £( 


-z 


( 1 . 2 . 11 ) 


k=0 


which is convergent in the circle \z\ < 1 in the complex plane. In that region, the sum is 

1 


P{z) = 


1 + z 


( 1 . 2 . 12 ) 


We can analytically extend the function P(z) to the entire complex plane, apart from z = —1. 
The starting series (1.2.8) is obtained when z = 1, outside the convergence domain, but with 
analytical continuation we can assign to the sum the value P = P(l) = 1/2. 


1.2.3 Borel summation 


Once that we have a divergent series, there are several ways to assign a value to its sum. In the 
following, we will concentrate on one method: the so-called Borel method. 

Let us define the Borel transform of a generic series, like that in equation (1.2.1), as 

OO 

fc =0 

In the Borel method, we define the sum of the series as 


S = 


dfe-‘B[P](t) 


(1.2.14) 


If the series S were convergent, we could exchange the sum with the integral and integrate term 
by term. By doing this, we would obtain that the sum of the series would be exactly S. 

We say that a series is Borel-summable B-summable if 
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• its Borel transform converges Vt 

• B[5](t) is defined on 0 < t < oo 

• the integral converges. 

We can generalize the Borel method and define higher-order Borel transform: 

/•oo roo roo 

Sn= dti dt2--- ... ,tn), (1.2.15) 

Jo Jo Jo 

where 

OO 

Bn[S]{ti,t2,. . ■ ,tn) = '^7^iJit2---tnt (1.2.16) 

k=0 ^ '' 

is the n-Borel transform. It can be shown that if a series is B,i-summable it is also B^-summable, 
Mk > n. 


1.2.4 Examples of Borel summation 

Once that we have defined the procedure of Borel summation, it is instructive to try and apply 
this method to some known divergent series. As we will see, most of the series that we are about 
to study have a physical meaning and correspond to a specific source of divergence. 


Alternating series 

The hrst example is just a simple exercise that confirms the results obtained for the previously 
cited series Eq. (1.2.8). The Borel transform of the series is 




k=0 


k\ 


The series converge in the whole complex plane, and the Borel sum is 


P = 


fOO 1 

dt 


(1.2.17) 


(1.2.18) 


like we previously found. 


Factorial divergence, alternating sign 

Consider the following series which, as we shall see, represents the divergent behaviour induced 
by the infrared renormalons: 

OO 

i?iR = J](-1)''A:!. (1.2.19) 

k=0 

Its Borel transform reads 

OO - 

B[/;ir](«) = = —, (1.2,20) 

k=0 

where the convergence radius of B„ for n = 1 is |t| < 1, while it can be shown that the n > 1 
Borel transforms have infinite convergence radius. The Borel sum is given by 

~ r°° 1 

Rik= / dte-' —. (1.2.21) 

Jo i +1 
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The result is the same at all orders 

Rj-^ = _eEi(-l) = 0.596347. (1.2.22) 

The function Ei(x) is called Exponential integral and is a special function defined on the complex 
plane as 

/ + 00 -t 

dt- (1.2.23) 

-X ^ 

The definition above can be used for positive values of x, but the integral has to be understood 
in terms of the Cauchy principal value due to the singularity of the integrand at zero. 


Factorial divergence, same sign 

Consider now the divergent series, similar to that of the previous example and physically related 
to the behaviour of ultraviolet renormalons; 


R^y = '^k\. (1.2.24) 

A:=0 

The corresponding Borel transform is 

B[7?uv](i) = (1.2.25) 

Once again, the first order Borel transform has convergence radius |t| < 1, and higher order 
transforms converge everywhere. However, for this series the Borel inversion intergal does not 
converge because of a pole {t = 1) in the integration path for n = 1 and because of the bad 
behaviour at t —)■ oo for the higher-order Borel transforms. However, the first order Borel integral 

Ruy= / dte-^- -- (1.2.26) 

Jo i — t 

can still have a meaning if we deform the integration contour in the complex t-plane and avoid 
the pole. But, for this reason, the result has an ambiguity, given by the two possible way of 
avoiding the pole (above or below the positive real axis). The result is 

iiuv = “ [Ei(l) ± ivr] . (1.2.27) 

e 

However, if we tried to use any of the n > 1 methods, we would end up with a different result. 
In this case, the only way to determine which result is correct is to consider the power series 

OO 

R{z) = ^{-z)^ k\. (1.2.28) 

fc =0 

1.3 Known sources of divergence 

We have seen that perturbative series in QCD are always divergent. We have also introduced 
some tools that can help us deal with divergent series. What is left to see is how the divergence 
arises in computations. There are three known sources of divergence in QCD that we are going 
to explore in this thesis: 
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• Landau pole divergence 

• Renormalons 

• Instantons 

The origin of the Landau pole divergence resides in soft gluon resummation. The basic idea is 
that soft gluon radiation has the effect of replacing the hard scale by a softer scale that is related 
to the process of radiation. For example, let us assume that a physical process is characterized 
by the hard scale and a generic scaling variable 0 < 2 < 1. What happens is that near the 
z = 1 region, close to threshold, the resummation of large logs of 1 — 2 ; replaces as{Q'^) with 
as(Q^(l — z)). Since ag has a rescaled argument, it becomes too large to be treated with the 
perturbative approach. The physical interpretation of this fact is that when 2 ; —)• 1, the center- 
of-mass energy is just sufficient to produce the given final state, so in this limit the process 
becomes elastic. 

In practice, at some low scale A (the position of the Landau pole) the strong coupling 
explodes, so when 

A2 

z = ZL^l- (1.3.1) 

resummed results become meaningless. The scale A is usually taken to be Aqcd- 

It is usually easier to consider resummed results in terms of the variable N which is Mellin 
conjugate to 2 ;. The 2 ; —)• 1 region corresponds to —>■ 00 . If the Ai-space resummed result is 
expanded perturbatively in powers of as{Q‘^) and then Mellin transformed back to 2 ; order by 
order, one ends up with a divergent series, and the source of this divergence is the presence of 
the Landau pole. 


Renormalons will be discussed in Chapter Firstly identified as inhnite chains of bubble 
diagrams, the renormalons are singularities in the Borel complex plane. They are related to 
small and large momentum behaviour of perturbative corrections, and divided respectively into 
Infrared (IR) renormalons and Ultraviolet (UV) renormalons. The divergent behaviour induced 
by renormalons strongly depends on the position of the poles. We will see that the Borel 
transform Eq. (1.2.13) of a perturbative series that has the renormalon problem presents poles 
at t = m/3o, with m an integer and /3o the first term of the /3 function. In particular, the closer 
the pole is to 0, the sooner (in terms of perturbative orders) the divergence occurs. There will 
then be a leading renormalon (m = 1) that will dictate the divergent behaviour of the series. 

Furthermore, we will see that renormalons are connected with the concept of power correc¬ 
tions. The fact that there are poles in the integration path of the Borel integral Eq (3.1.12), 
leads to an ambiguity in the results. This ambiguity can be interpreted in terms of the differ¬ 
ence between the asymptotic value of the series and the real sum. We will see how these power 
correction vary with the scale and how this scaling behaviour is related to the position of the 
leading renormalon. 


Instantons are field configurations fulfilling the classical equations of motion in Euclidean 
spacetime, which can be interpreted as a tunneling effect between different topological vacua. In 
quantum mechanics there are cases, e.g., the double-well potential, in which the ambiguity intro¬ 
duced by the Borel summation can be cured by a procedure called the Bogomolny-Zinn-Justin 
(BZJ) prescription. Here the perturbative ambiguity cancels against a non-perturbative contri¬ 
bution from instanton-anti-instanton events. The sum of the perturbative and non-perturbative 
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semiclassical expansions in quantum mechanics apparently produces ambiguity free (and accu¬ 
rate) results. 

One may ask if this idea can work in field theory, and in particular in QCD. The answer is 
that it does not work for gauge theories on due to the above-mentioned IR renormalon prob¬ 
lem. Reference argues that this method does work on x in a gauge theory continuously 
connected to one on but the generalization to QCD is an open problem. 



Chapter 2 


Divergence in resummed series 


In this Chapter we study one of the known sources of divergence in perturbative QCD, the Landau 
pole divergence. First, we will review the basics of soft-gluon resummation. Our discussion will 
mainly follow Refs. 


We have already anticipated in Section O that, because of the presence of the Landau pole, 
the inverse Mellin transform of the resummed result in the soft-gluon limit does not exist. If we 
expand the resummed cross section in powers of q;<j, each of the terms has a finite inverse Mellin 
transform, but the correspondent perturbative series in physical space is divergent. We will 
describe a prescription to obtain the resummed result in 2 space, the Borel prescription [^. This 
is not the only way to deal with the Landau pole divergence: another widely used prescription, 
which we will not cover here, is the minimal prescription [^. 

Finally, we will show that an approximation can be constructed using the expansion of the 
resummed cross section in the soft limit. For processes like Higgs production in gluon fusion, this 
approximation succeeds in predicting fixed order calculations. The divergent behaviour of the 
expansion of the resummed cross section is related to that of the exact fixed order calculation. 


2.1 Soft-gluon resummation 

The generic coefficient function C{z,as) suffers from kinematical enhancements due to gluon 
emissions when 2 ; —)• 1. The divergences that appear when calculating virtual contributions are 
exactly cancelled by those soft gluon emissions. In fact, a generic physical cross section needs 
to be inclusive over arbitrarily soft particles in the final state, since any detector has a finite 
energy resolution. Nevertheless, soft-gluon effects can still be large in some kinematic regions. 
Because of that, calculations to all orders of perturbation theory are necessary to achieve reliable 
predictions. In this framework, resummation is basically an all-order summation of certain 
classes of logarithms. 

To see this in practice, consider as an example the case of a quark parton line emitting n 
gluons, as in Fig. 2.1, The energy fraction of each emitted gluon to the quark energy is 1 — Zi. 
Therefore, after n emissions the energy of the quark will be 2 : = 2 : 12:2 ••• 2 ;^ times its initial 
energy. When we integrate over the phase space of the emitted gluons, we have a kinematic 
enhancement for each gluon. It can be shown that all these terms convert into the sequence 


a. 


log"(l - 

1-z 


0 < k < an — 1 


( 2 . 1 . 1 ) 
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1 - Zi 1 - Z2 1 - Z3 1- Zn 

Figure 2.1: Emission of n gluons from a quark parton line. The quantity zt represents the energy 
fraction for each line. Figure taken from Ref. 1^. 


where a = 1 for deep inelastic scattering and a = 2 for Drell-Yan and Higgs production (we will 
consider this last case). 

It is clear now why we need resummation: there will always be a certain kinematical region 
where z is such that 

a,log2(l-z) ~ 1. (2.1.2) 

In this case, all terms in the perturbative series are of the same order, and any truncation would 
neglect a huge part of the result. 


2.1.1 Resummation in Mellin space 

The concept of resummation has been known for a long time |9 - 11 . This section only serves as 
a reminder of the most important concepts. 

Whenever we try to compute a cross section, we need of course a matrix element and an 
integral over the phase space. In the case of n gluon emission from a parton line, the calculation 
of the matrix element can be performed in the eikonal approximation, according to which the 
matrix element dJln factorizes as 


1 

Tln{zi, . . . , Zn) ""zz — TT (2;^), 

n\ 


(2.1.3) 


Z=1 


where Dili is the matrix element for the single emission. Unluckily, the phase space in physical 
space is not factorized, because of the Dirac 6 that expresses the conservation of momentum 


dzi dZ 2 ■ ■ • dZnS{z - Z1Z2 ■ ■ ■ Zn). 


(2.1.4) 


However, in Mellin space even the phase space factorizes 

f — Z^ S{Z - Zi---Zn) = Z^~^ ■■■ Zn~^. (2.1.5) 

Jo Z 

We remind that the Mellin transform of a function f(z) is defined as 

fiN) = M[f]iN) = C dz z^-^fiz). (2.1.6) 

Jo 

The Mellin transform can be seen as a Laplace transform where a change of variables has been 
performed. The inverse Mellin transform is 


f{z) 


1 

27rz 


rc-\- 

J c—i' 


C+200 


dN z-^fiN), 


(2.1.7) 
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where c is greater than the real part of the rightmost singularity. It is easy to show that the 
Mellin transform of a convolution is the product of the Mellin tranforms 


if^g){N) = dx X 
Jo 


N-l 


dy 


r-l 


—f{y)g {-]= I dx X 
y \yj Jo 

rl /■! 


N-l 


= / dy y^^-^f{y) / dzz^^-^g{z) = f{N)g{N). 
Jo Jo 


dy dzf{y)g{z)5{x - yz) 

Jo 

( 2 . 1 . 8 ) 


Back to the emission of multiple gluons, it can be shown that the threshold region z ~ 1 
corresponds to the region of large Mellin moment N. In particular, the sequence of contributions 
Eq. (2.1.1) converts into the tower 


a” log —, 0<k<2n. 

Therefore, the coefficient function in iV-space at order a” is given by 

ci")(iv) i [c« (jv) 

where is the Mellin transform of soft terms up to order If, for example, we only 


(2.1.9) 


( 2 . 1 . 10 ) 


'soft 

consider the leading term, the coefficient function becomes 


C. 


( 1 ) 

soft 


(iV) = [ dzz 

Jo 


N-l 


4^1 




( 2 . 1 . 11 ) 


where Ai = CpjTi for the Drell-Yan case and A\ = CaJti for the Higgs case. 

Now it comes naturally that the soft terms can be resummed explicitly, and the result is the 
exponential: 




n=0 


- soft 


= exp 


“.cW(iv) 


( 2 . 1 . 12 ) 


Since we have only considered the first order in as, we say that this result is valid only at 
leading-logarithmic (LL) accuracy. In fact, only the highest power k = 2n in the tower of logs of 
Eq. (2.1.9) is resummed. Eurthermore, this expression does not take into account the running 
of q;<j. It can be proven that the most general expression for the iV-space resummed coefficient 
function is 


- / 

= go{as)expS ^ 

where S is called Sudakov form factor and is defined as 


m2 


fi 




1 


^ — z Jm^ 


_/Vf2(i_2)2 2 

-^2H {as{y^)) + D (a,([l - zfM^)) 




(2.1.13) 


(2.1.14) 


We have introduced the functions go{as), A{as) and D{as). They are always represented as 
power series in Og, with 5o(«s) = 1 + 0{as). D{as) and go(as) are process-dependent, while 
H(as) is process independent: it is the coefficient of the soft singularities in the Altarelli-Parisi 
splitting function. The resummed coefficient function can be written in a form where the Sudakov 
exponent is organized in powers of logs and as ■ 


m 2 ) = go{as) expS {aL,a) 


( 2 . 1 . 15 ) 
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where qq can be written as 


a = 2/3oa, L = log—, 


(2.1.16) 


9o{as) = 1 + '^90jai, 

and S has the following logarithmic expansion 

S(aL, a) = —gi{aL) + g 2 (aL) + ag^iaL) + a'^gJaL) + 
a 


(2.1.17) 


(2.1.18) 


The functions gi can be obtained performing the integrals Eq. (2.1.15), and they are determined 


(as shown, for example, in Refs. 00 ) by a small number of coefficients of the expansion of 
A(as) and D{as), and are of order gi{aL) = 0{ag), gi{aL) = 0(as) for f > 1. 

One finds that the N^LO anomalous dimension is necessary in order to obtain < 7 ^+ 1 , which 
enters the result at N^’LL accuracy. Conversely, most of the times the process-dependent func¬ 
tions are determined by matching the expansion of the resummed result with a fixed order 
computation. Finally, predictions for phenomenology at N^LO-|-N*'LL accuracy are obtained by 
combining the the fixed-order computation with the resummed coefficient function expanded in 
powers of a^. In this step it is important to remember to subtract the double-counting terms: 


^NPLO 

^N'^LL 


(iV,a,) = + 

j=0 


Ops 


i=o 


N'^LL'- 


jj fives 

® '-^N''LL 


(fV, as 


dai 


(2.1.19) 


j 03=0 


2.2 The Landau pole divergence 

So far we have seen how soft gluon resummation is performed in Mellin space. As we anticipated 
in Section |1.3[ the resummed cross section in Mellin space has a branch cut for large values of N 
(soft region). This means that if we expand in ag the resummed cross section and then compute 
the inverse Mellin transform term by term, then each inverse transform exists, but the resulting 
series is divergent in z space. Let us see how this divergence arises and why it is related to the 
Landau pole. Consider a generic observable in physical space cr{Q‘^, z) and its V-space transform 


a{Q‘^,N)= / dzz’^ ^a{Q‘^,z). 

Jo 


( 2 . 2 . 1 ) 


Since the resummation has the form of an exponentiation, it easier to work in terms of the 
physical anomalous dimension, which is defined to be 


l{ots{Q‘^),N) = 


d\D.a{Q‘^, N) 
d\nQ‘^ 


( 2 . 2 . 2 ) 


Let us consider for now the structure functions for DIS. Similar observations can be made for 
most of the other processes that are relevant at the LHC (e.g. DY, Higgs production). The 
resummed expression of the physical anomalous dimension ^{as{Q‘^).iN) has the form [9- 


11 


rN 


dn 


7(a.(g2),iV)= / —Y.9k»siQVn) + 0{N^), 

./i n 


(2.2.3) 


k=l 
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where gk are constants which be determined either from first principles [9,10 or by comparison 
with the fixed-order calculations. Here, the sum over k represents the sum over successive 
orders in logarithmic accuracy. The resummed expression of 'y{as{Q‘^), N) at N^LL can be used 
to compute the resummed cross section in Mellin space at N^LL, but then one needs to compute 
the inverse Mellin transform of this quantity to obtain the resummed cross section in physical 
space. 

The main effect of resummation, as can be seen explicitly in Equation (2.2.3), is to rescale 
the argument of the strong coupling to a softer scale as{Q‘^/N). As we are going to show, this 
replacement corresponds to the introduction of a branch cut in the positive real A^-axis, which 
makes it impossible to perform the inverse Mellin integral. For this reason, 7 (as(( 5 ^), A^) cannot 
be transformed back in physical space. Let’s see why and where the branch cut occurs in the 
simpler case of LL accuracy: 7 has the form 


r-iV 


dn 


'yLL{as{Q'^),N) = gi I —as{Q'‘‘/n) = -^ln{l + /3oas{Q'‘‘)ln^ 


1 


N 


Jl Po 

To be consistent, remember that we need to use the the leading-log expression of 

as{Q‘^) 


as{g ) = 


1 


+ PoOisiQ'^) log ^ 


(2.2.4) 


(2.2.5) 


If we take Q^jn as the argument of as, wee see that the denominator of Eq. (2.2.5) becomes 


1 - (3oas{Q^)logn. 


( 2 . 2 . 6 ) 


The singularity starts when this denominator vanishes, namely 'yhhioisiQ'^), Ai) has a branch cut 
on the real positive axis for values of N that satisfy 


N > NXj = ^ 


(2.2.7) 


From the presence of the branch cut follows that the inverse Mellin transform of the physical 
anomalous dimension does not exist. However, one may consider the inverse Mellin transform 
of each term of the expansion of 7 ll(q:s(<3^)) N) in powers of as{Q‘^)- This would be 


( 1 ^ I 1 T 


rN-\-ioo 


dNx-^ In^ ^ 


N—ioo 


1 _ 

iV’ 


N > 0 . 


( 2 . 2 . 8 ) 

In this case, each term of the series is a well-defined inverse Mellin transform, but the series as a 
whole is divergent, so the limit AT —)• 00 is not defined. If the series were convergent, one would 
be allowed to interchange the sum and the integral in eq. ( 2 . 2 . 8 ), but then the sum 


(2.2.9) 


would only be convergent when 


Poas{Q^)ln 


1 

N 


< 1 . 


( 2 . 2 . 10 ) 


There are several ways to deal with the Landau pole divergence. Here we will discuss the 
so-called Borel prescription. Another way is what is usually referred to as the minimal prescrip¬ 
tion [^, which we will not cover in this thesis. 
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2.3 Borel prescription 


In order to understand the Borel prescription method, let us work with the divergent perturbative 
series eq. (2.2.8). Following step by step Ref. [^, we compute the Mellin inversion integral: 


1 pN-\-ioo jk 1 pN-^ioo 

-: / dNx-^ — : dNx-^N-^ 

27ri JN—ioo 27ri jN—ioo 


d^ 1 


d^k 


In^-i - 


X 

+ 


ri=0 

+ (5(1 - x) 


ri=0 


where we have used the identity 

/■i 


dx X 


N-l 


10 


In^ 


-1 


= r(r/)(iV-^-l). 


Hence we obtain 
PLL{asiQ^),x) = 


(2.3.1) 


Po 


E 


k=0 


k + l 
+(5(1 - x) > = 


E 

I n=0 


k + l \ 

f d^ 

1 ^ 

^k+l-n 

- 


n J 

\dq^ 

nv)j 

dqk+^-n- 

X 

+ 


ri=0 


9l \~PoOts{Q‘^)]^~^^ 


/3o 


k + l 


k=0 

+(5(1 - x)| = 

9l q „ ,'^2',in+l / A("')(1) 1 

-2^[-/3oa.(Q )] <1^^ ^2^ 


X n=l ^ ^ ^ ^ 




/3o 


n=0 


1 


X k=n ' ' 


+ —-(5(1 -x) ^ = 
n + I 


K 


+0{a 


^[-/3oa,(Q2lni)]-+i 

In - X 


n + 1 ^ ^ 




(2.3.2) 


Here we have defined A{z) = l/r(z), and we have used the identity A^^\0) = kA^^~^\l). In 
the limit R —?• oo the terms of order can be neglected, but we have already seen that the 

series is divergent. 


In the large x limit and up to leading logarithmic accuracy we can rewrite eq. (2.3.2) as 

' R{as{Q‘^),x 


Phh{as{Q‘^),x) = ^ 
Po 


1 — X 


(2.3.3) 
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where 


K 

RiasiQ^),x) = hm (2.3.4) 

K^oo ^^ 
n=0 


The idea of Ref. is to sum the divergent series using the Borel method that we introduced 
in Section [l. 2.3 Namely, we take the Borel transform of the divergent series (2.3.4) with respect 


to f5oasiQ {I - x)): 


R{w, x) = -Y^ -^ {-wy = - 


1 


j=0 


T{l-w)' 


(2.3.5) 


To obtain the sum of the series, we have to compute the Borel integral 


RB{asiQ‘^),x) = - 


r+oo 

Jo 


dw e Po°‘a(Q^(i-x)) 


1 


10 r(l-u;)' 

Unfortunately, the integrand diverges as tc —)• oo. In particular, the reflection formula 


1 


r(l — w) TT 


= —T{w) sin(7rt(;) 


(2.3.6) 


(2.3.7) 


implies that A(1 — w) has a factorial oscillating behaviour when w —)> oo. 

We already saw in Section 1.2.4| how to deal with singularities that are along the path of 
integration in the Borel inversion integral. In our case, the singularity is at tc —>• oo, hence 


we simply introduce an upper cutoff C to the integral. The divergent result eq. (2.3.6) is then 
replaced by 


rC 


1 


RB{asiQ),x,C) = - dw 

Jo r(l-ri;) 




(2.3.8) 


which is convergent for all C and well defined for all x. Indeed, we can expand the integrand 


using eq. (2.3.5) and integrate term by term. The result is a convergent series, but, if (7 —?• oo. 


then lim /fc = 1 and we obtain once again the original divergent series. 

C^oo 


2.4 Soft-gluon approximation 

We have seen how the expansion of the resummed cross section in z space is a divergent series. In 
some cases, at low perturbative orders the terms obtained by expansion of the resummed result 
turn out to provide a good approximation to the full result, even far from threshold. In the case 


of Higgs production, this has been known for long. Moreover, in Ref. 12 it was shown how the 


information from soft resummation can be combined with that from high-energy resummation 
to construct an approximation 


C*approx(.^) Ols) — C*soft(N, ds) -|- C'h.e.(N, O^). 


(2.4.1) 


For the purpose of our thesis, we will only be interested in the soft part of this approximation. 
To see how we can use the resummed result to predict the full fixed-order computation, let us 
follow Ref. [12] . We first compare the resummed coefficient function in N space with the exact 
expression CJdo{z), which is given by 

= 4Ag(z)'Di{z) +d6{l - z) - 2Ag{z)^^^ +TZgg{z), 


( 2 . 4 . 2 ) 
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where 


, , f ln*^(l — z) 


(2.4.3) 


and 


Ag{z) = 


Ca1-2z + 3z^ - 2z^ + 


(2.4.4) 


vr z 

Here, d and TZgg{z) are functions of the dimensionless ratio niH/mt. The function TZgg{z) is 
regular in the limit iV —>• oo, and therefore we won’t need to worry about it. 

The counterpart of this comparison is the expansion of Eq. (2.1.15) to 0{as) at NLL loga¬ 
rithmic accuracy: 


C*res {N, as) = l + a.C«(Af) + 0(0^), 

(iV) = ffi ,2 ln2 N + 52,1 In + 5o,i, 


with 


51,2 = 


2Ca 


IT 


4Ca 

52,1 =-7e- 

vr 


Now, if we compute the inverse Mellin transform of Eq. (2.4.6), we get 

Crei (z, as) = 5o,i<5(1 - z)+ 25p 2T>5°®(2:) -h {2-fEgi,2 - 52,i) T>o°®(z), 


= 50,1(5(1 - z) + 


vr 


1 (z), 


where 




'in^ inf 
Inl 


(2.4.5) 

(2.4.6) 

(2.4.7) 


(2.4.8) 


(2.4.9) 


If we now compare (2.4.9) to the soft contribution to the exact coefficient function Eq. (2.4.3), 
we notice that they are different. The reason of this difference can be explained by noticing 
that singular terms as z ^ 1 appear in the integral of the real emission diagrams over the gluon 
transverse momentum. 

The approximation worked out by Ref. 12 is different from Eq.( 2.4.8) by subleading terms. 
The validity of the approximation can be tested by comparing the soft approximation to the 
full result. Results as functions of the Mellin moment N are shown in Fig. 2.2 for the hrst two 
perturbative orders of the coefficient function. 
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Higgs partonic gg. Soft approximations. Order Og Meiiin transform 



2 

Higgs partonic gg. Soft approximations. Order Og Meiiin transform 



Figure 2.2: Plot of the partonic coefficient functions and C^‘^\N) for mn = 125 GeV. 

The black line represents the exact result, while the soft approximation is presented in various 
forms; the two preferred approximation of Ref. |12| are called softi and soft 2 , the so-called N- 
soft approximation (the one based on as in Eq. (2.4.9)), a collinear-improved A^-soft 

approximation worked out by Ref. [I^ and the soft-0 approximation which we didn’t discuss. 
Figure taken by 
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Chapter 3 

Renormalons 


Another source of divergence in QFT is known as renormalons. The ultimate origin of the 
divergence caused by renormalons can be found in the large momentum and small momentum 
behaviour of certain classes of corrections that can be inspected at all orders in perturbation 
theory. Giving a definition to the concept of renormalon is not easy, so we will start with a 
classical example, the bubble diagram chain, in order to introduce the subject. A complete and 
exhaustive review on renormalons can be found in [14| . After introducing the issue, we will show 
that there is a deep connection between the renormalons and the complex plane where the Borel 
transform (that we saw in Section 1.2.3) is defined. Finally, we will explore the Borel plane 
and list all the known singularities that need to be taken into account when performing Borel 
summation. 


3.1 The bubble diagram chain 

We will now introduce the idea of renormalon divergence Ref. [14] . We will start with the 
computation of the bubble diagram chain and then we will give it an interpretation in terms of 
renormalons. 

In particular, we will consider the case of the correlation function of two currents of massless 
quarks 




d^xe (0|r (j^(x)>(0))|0) = {q^,qu - 


(3.1.i: 


where, as usual, = —cp'. We will consider a certain class of corrections to the Adler function, 
which is defined as 

dn(Q2) 


D{Q^)= 47r" 


dQ^ 


(3.1.2) 


The class of corrections that we would like to compute is represented in Fig. 3.1 Keep in 


mind that this is just an illustrative example: they don’t represent the only contribution to the 
renormalon divergence. However, historically, this is how renormalons were first introduced in 
the original works by Refs. [4,15,^. To compute the Adler function we need to insert multiple 
times the renormalized fermion loop; 


/3o/«s [ln(-fc^//x^) + C] , 


(3.1.3) 
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Figure 3.1: The set of bubble diagrams for the Adler function consists of all diagrams with any 
number of fermion loops inserted into a single gluon line. Figure taken from Ref. 
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where C is constant that depend on the renormalization sheme [C = —5/3 in the MS scheme) 
and /3o/ is the fermion contribution to /3o 

/3o, = -g. (3.1.4) 

For the time being, we are only considering fermion bubbles, so we only need the fermion 
contribution to the (3 function. We will see what happens when we include gluon and ghost 
terms in Section [3.1.11 

When we want to include n —)• oo fermion loops, we integrate over the loop momentum of 
the big external fermion loop and over the angles of k. Let us define = —k‘^jQ‘^, we get 


“ f Jt2 


n=0 




/3o/asln k 


,Q^e 




(3.1.5) 


The function F is rather complicated in its exact form, but we can just consider the n ^ 
1 approximation, due to large logarithmic enhancements, the dominant contributions to the 
integral come from very large or very low values of /c (A: ;g> Q and k Q). For this reason, we 
are only interested in the small-A: and large-fc behaviour of F: 

F{P) = —k^+ 0{k^ In k'^) k^O, (3.1.6) 

27r 


It is worth while to notice that we expect a power-like approach to zero in both cases since the 
Adler function is UV and IR finite. If we split the integral (3.1.5) at = fj?j and we 
insert (3.1.6) for the small-fc^ interval and (3.1.7) for the large-A: ^ interval, we obtain 


n =0 


/g2 5\-2 /pofY 1 g2 


n! 


11 

n+ — 

6 


(3.1.8) 


where the first term comes from the small k region and is accurate up to relative corrections of 
order n (2/3)”, while the second term comes from the large k region and is accurate up to (1/2)”. 
This is the reason why the factorial divergences of the two series are called respectively infrared 
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(IR) renormalon and ultraviolet (UV) renormalon. The etymology of the word “renormalon” 
goes back to Ref. [^. There, the word renormalon was chosen as an analogy to the only other 
known (at that time) source of divergence: the instanton divergence. This divergent behaviour 
was then new and typical of renormalizable field theories. 


Equation (3.1.8) exhibits the same factorial divergence that we encountered in the examples 


of Section 11.2.41 The first order Borel transform reads 


Blow = 


_5 
ye 3 


-2 


1 


2 - Poft 


Cf _5 
+ X 2 ® ^ 

dvr 


5 

+ X- 


(l + /3o/t)2 6(l + /3o/t) 


(hrst IR renormalon) 


(first UV renormalon), (3.1.9) 


The corresponding singularities in the Borel plane lie at t = 2//3o/ (IR renormalon) and t = 


— l/Pof (UV renormalon). Eq. (3.1.9) only gives us the singularities that are close to the origin. 
It can be shown that the exact Borel transform of the set of diagrams of Fig. |3.1| is made of an 
infinite sequence of IR (UV) renormalon poles at negative (positive) integer multiples of m/3o/t, 
apart from m = 1. Therefore, we define the term renormalon more generally as a singularity of 
the Borel transform, which is ultimately related to the large or the small momentum behaviour 
of the loop. As we have specified above, the set of bubble graphs only provides some of such 
singularities, but not all of them. 

Now that we have introduced the subject of renormalon divergence with an explicit exam¬ 
ple, let us make a few observations, which will be useful for the application of renormalons to 
determine theoretical uncertainties in QCD. 


3.1.1 Adding gluon and ghost terms 


So far we have only considered one loop corrections due to fermion bubble graphs, and therefore 
we have consistently used the fermion contribution /3o/ to the /3-function. The next step is 
to add the gluon and ghost bubbles, but that introduces a complication, since the result is 
gauge-dependent. 

What happens is that the effect of substituting /3o/ with /3o flips the location of the renormalon 
singularities. The fermion contribution /3o/ is negative, but /3o is positive, and therefore UV 
renormalons are now located in the negative real axis, while IR renormalons are in the positive 
real axis of the Borel plane. Moreover, UV renormalons give origin now to a sign-alternating 
divergence, while IR renormalons introduce an ambiguity in the Borel integral as we saw in 
Section 1.2.4 To remove this ambiguity, one needs to add non-perturbative corrections. The 


same situation occurs in QED with UV renormalons. In both cases, as one could easily expect, 
non-perturbative corrections are needed when the coupling becomes large (infrared in QCD and 
ultraviolet in QED). 

At this stage, there is no evident reason to suppose that we can extrapolate to the full non- 
abelian /3o. However, it is shown by that the substitution of /3o/ by /3o can be fully justihed. 
For further simplihcation, since the poles of the Borel transform are located at multiples of Pot, 
we will use the dehnition 

u = Pot. (3.1.10) 


This way, IR renormalons are located at u = n, while UV renormalons are located at n = —n, 
where n is a positive integer 7 ^ 0 . 
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3.1.2 Power corrections 


Let us now consider a generic perturbative series where we have already collected an overall as 
factor 


R = '^rka 


k+1 

s 


k=0 


This way the Borel integral has a slightly different form than in section 1.2.3 

roo 

R= dte-*/“*B[i?](t). 

Jo 

If the Borel transform B[R]{t) presents the typical IR renormalon behaviour, 

K 


B[R]{t) = 


1 — aPot ’ 


(3.1.11) 


(3.1.12) 


(3.1.13) 


where 1/a is a positive integer and the radius of convergence is |t| < l/|a/3o|j then we have 
an ambiguity in the determination of the integral. R acquires an imaginary part, which is the 
residue at the pole of the Borel integral: 


lm(i?) = TvrRes ( ) = ±7r—— exp (--- ) . (3.1.14) 

\al3oJ a/3o V ajjoas J 

Now we might be interested in the scaling behaviour of this ambiguity, which in some sense 
represents the difference between the exact value of the series and the asymptotic value given 
by the Borel integral. To do so, let us recall the definition of the running coupling constant in 
terms of Aqcd: 

asii?) = -(3.1.15) 

/3o log -72- 

^'-QCD 


If the Borel transform has a singularity at u = 1 /a, that yields an ambiguity in the definition of 
the Adler function which scales as 


6D{Q^) oc exp ( — 


a/3oas{Q'^) 



( Aqcd \ “ 

V Q ) 


(3.1.16) 


Recall that the first IR renormalon pole in (3.1.9) corresponds to a 
leading power correction proportional to 


1/2, and therefore to a 




(3.1.17) 


3.2 The Borel plane 


We have seen how a renormalon is nothing but a pole on 
integral. For what concerns the Adler function, the Borel 
course, there might be a whole new set of singularities that 
known singularities we distinguish three sets: 


the real axis of the Borel inverse 
plane is shown in Figure 3.2 Of 


we don’t know, but among all the 
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t 


Instanton-anti-instanton 
singularities at t = 47r, Stt, ... 



U\" renornialons 

at f = —7ti/do, Jii = 2,3,... 


IR renornialons 

at / = m/fio, '» = 2 


Figure 3.2: Singularities in the Borel plane of the Adler function in QCD. The singular points 
are shown, but not the cuts attached to each of them. Figure taken from Ref. 


14 


UV renormalons are poles located in the negative real axis at t = —mjl3Q, with positive 
integer m. The leading term of the power correction associated with UV renormalon should be, 
following (3.1.9), of order using (3.1.16). Actually, it can be shown that its precise 

form is 


SDuy 


g2A2 




X logarithms. 


(3.2.1) 


However, the Borel inverse integral only involves positive values of u, so UV renormalons produce 
no ambiguity in the Borel integral and a sign-alternating factorial divergence. Therefore, no extra 
terms should be added to the perturbative expansion. UV renormalons depend on the theory, 
but are process independent. 

IR renormalons can be found in the positive real axis, at t = m//3o, with m = 2, 3,.... The 
leading term associated with the first IR renormalon is of order (A/Q)^. Unlike UV renormalons, 
IR renormalons are process-dependent. 


3.2.1 The instanton divergence 


In Figure 3.2 we can observe a third set of singularities, which leads to a factorial divergence 
in the perturbative series: the so-called instantons [17] . Instantons are classical solutions to 
equations of motion with a finite, non-zero action. Ref. 18 shows how configurations of n 


instantons and n anti-instantons give rise to poles in the positive real axis in the Borel plane 
located at t = Ann. Instanton divergence can be related to the factorially increasing number of 
diagrams in perturbation theory. Thanks to the semi-classical origin of instantons, the residues 
at the poles in the Borel plane can be calculated. This calculation is worked out by 19 for the 


Adler function. However, in QCD instanton poles are located at large t, where t is the variable of 


the Borel transform. We shall see in Section 5.4 that this means that their effect is by all means 


negligible for what concerns large orders in perturbation theory in QCD. Moreover, recalling 
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Eq. (3.1.16), it is easy to see that they do not represent a dominant source of power corrections. 


For a complete review on the subject of instanton, see 20 





Chapter 4 

Theoretical uncertainties 


Whenever performing an experiment, such as at the LHC, we compare measurements to the¬ 
oretical calculations and try to find out if they match or not. However, in QFT the theory is 
perturbative, which means that theoretical predictions are perturbative series of which only the 
first terms are known. In this case, a full control of the uncertainty of these predictions becomes 
of paramount importance, as both the experiment and the theory need to be provided with a 
degree of uncertainty in order to determine their agreement. In QCD the issue of theoretical 
accuracy is pressing, due to the large size of the coupling ag and therefore its slow perturbative 
convergence. 

There must exist, therefore, a definition of theoretical uncertainty on any calculation of 
any observable at the LHC. Following Ref. 21 , we first make a distinction between parametric 


uncertainties (PU), related to the value of input parameters, and actual theoretical uncertainties 
(THU), related to our lack of knowledge about higher orders in perturbation theory. There is no 
way of eliminating completely parametric uncertainties, even though they can be reduced when 
more precise experiments produce improved results. Theoretical uncertainties, however, might 
be ideally eliminated if an all-order computation were available. 

Another difference between PU and THU is that PU are distributed according to a known 
(usually Gaussian) distribution while THU are arguably distributed according to a flat distri¬ 
bution, even if the statistical interpretation of THU is less clear. PU have been studied in detail 
during the past several years, and we will not talk about them in this thesis. Among the mod¬ 
els that give an estimate to the THU, we will consider three: the scale variation method, the 
Cacciari-Houdeau method and the David-Passarino method. 

The scale variation is the conventional and most widely used method to determine theoretical 
uncertainties. It is based on the following idea: in the full theory there should not be any scale 
dependence and order by order in perturbation theory we should be able to see the asymptotic 
limit. Therefore, variation of the scale (or scales) is a pragmatic way of understanding how far 
we are from controlling the theory. 

The Cacciari-Houdeau method, introduced in Ref. |^, is based on a Bayesian model that, 
given a certain distribution for the coefficient of the series, allows one to characterise a pertur¬ 
bative theoretical uncertainty in terms of a credibility interval for the remainder of the series. 

The David-Passarino method j^, instead, uses the concept of sequence transformation to 
improve the convergence of the series. Using sequence transformation, slowly convergent series 
can be transformed into series that have better numerical properties. Even if a sequence trans¬ 
formation hardly ever sums a series exactly, it usually predicts some of the unknown terms of 
the sequence. 
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In order to fix the notation, from now on we will define the partial sum of a perturbative 
cross section up to the k order as 


^ ^ CnOtg , 
n=l 


(4.0.1) 


and the remainder of the series as 


OO 

Afc = 'y ^ CnOtg . 
n=fc+l 


(4.0.2) 


4.1 Scale variation 

In QCD, there is no obvious optimal choice for the renormalization and factorization scale. In 
QED, we have a physical subtraction point, = 0, for photons with momentum transfer q, 
which is referred to as the Thomson limit. In the EW theory, once again, we have a physical 
subtraction point(s): the electromagnetic coupling is still fixed in the Thomson limit, while the 
weak mixing angle is related to the ratio of the W and Z boson masses. In other words, in 
the electromagnetic and EW theory our calculations do depend on hr, but this dependence 
disappears once the Lagrangian parameters are replaced by data. 

However, in QCD we have no analogue of Gi;’. This means that if s is the scale at which we 
study the process, our LO calculation will always contain logarithms like \n{s/HR)- One should 
hnd a scale where some data is available (a subtraction point), but this is not yet possible in 
QCD. So the question arises of what is the best choice for hr^ 

The general guideline 21 is to set hr exactly to the relevant scale s. Of course, this is 


straightforward in processes where only one scale is relevant, but in processes where multiple 
scales are relevant there will be additional logarithms of argument s/s'. In this case, the conven¬ 
tion is to choose the renormalization scale hr the factorization scale fiR process by process, 
in such a way to minimize the effect of the new corrections when going to the next perturbative 
order. This choice of scale is sometimes called dynamical scale and it is the closest thing to a 
subtraction point in QCD. 

The scale variation method consists in giving an estimate to the THU in terms of an interval 
centered around the dynamical scale, e.g. s/n < hr,f < ns. There is no specific choice for the 
value of n, but it should be fixed so as to include a plateau in the scale dependence. 

When combining the renormalization scale and factorization scale, once again we should 
try to minimize the effect of new corrections. The possible choices for this combination are a 
diagonal scan, a diagonal scan with anti-diagonal corrections or a two dimensional scan with 
1/n < hr/iJ-f < n. 

4.2 Cacciari-Houdeau model 


The main shortcoming of the scale variation is that it does not provide the degree of belief of 
the resulting uncertainty bands. In other words, it does not associate a numerical value to the 
probability that the uncertainty band contains the true sum of the series. Because of that, it 
becomes difficult to combine THU with the above mentioned PU. 

The Cacciari-Houdeau Method 22 provides theoretical uncertainties with a well defined credi¬ 


bility measure and computes explicitly the degree of belief of a given interval. 
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Of course, this degree of belief has nothing to do with the concept of frequentist probability, 
but it needs to be intended in terms of the Bayesian probability. While frequentist probability 
is linked to a large number of realizations of an experiment, Bayesian probability deals with the 
mathematical treatment of the “trueness” of a statement. The variables that allow us to treat 
the frequentist probability are called random variables, while we can call uneertain variables the 
ones that appear in a credibility distribution. 


4.2.1 Model overview 


Given a series like in Eq. (4.0.1), the model gives a form to the distribution of the coefficients 
Co, Cl,.... The three main hypotheses of the CH models are 


1. The residual density for the unknown coefficients Cn, given a known upper bound c, is a 
fiat distribution. 

2. The upper bound c is the only parameter that contains all the information. 

3. All values for the parameter c are equally probable. 


Once we have these three hypotheses, we have defined the credibility measure over the space of 
uncertain variables {cq, ci, ... }. 


Let us consider now a perturbative series whose first order is not n = 0 but n = 1. Using 
the three hypotheses, and defining = max(|ci|,..., |cfc|), the conditional density for the 
remainder can be obtained (see 22 for a complete derivation): 




E 


(y. o C-Y) 


n=k-\-l 


/(cfc+l, Cfc+2) • • • \Ch ■ ■ ■ ,Ck) dCfc+i dCfc +2 .... 


(4.2.1) 

It is difficult to treat this expression analytically, but we can make the approximation that the 
whole remainder of the series is comparable to the first term of the remainder itself: 


|Afc| ~ a|+^|cA:+l| . 

In this case the number of known coefficients is Uc = k + 1 — I, and we obtain 

1 


f{Ak\ci ,... ,Cfc) ~ 


nr_ 


1 


if |Afc| < a5+^C(A,) 


?T'c + 1 / 2a»^^Ci 


" [ (|A,|/(a5+l5(,)))-c+i 


if |Ai,| > 


C(k) 


(4.2.2) 


(4.2.3) 


If we know /(Afc|c;,..., Cfc), we can calculate any p-credible interval for A^, where p is the 
Bayesian probability that the true remainder of the series will be inside the interval: 


P= . ^ /(^fclQ, • •. ,Cfc)dAfc . 


(4.2.4) 


Finally, using a suitable approximation, it can be shown that 


4 '” = 


“S Hk) ric 


p 


if p < 


+ 1)(1 - P)] if P 


> 


ric+l 


ric+l 


(4.2.5) 
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where p = P/100 and P is a number between 0 and 100. 

Therefore, the CH model not only gives us the theoretical uncertainty interval for a generic 
perturbative series, but it also tells us what is the Bayesian conhdence level of that interval. 
However, the first hypothesis is that the coefficient of the series Cn should be of the same order 
of magnitude. This turns out to be not applicable to series like the Higgs cross section (see 
Table 1.1). The modified Cacciari-Houdeau approach, introduced in 23 , aims to solve this 
issue. 


4.2.2 The modified Cacciari-Houdeau approach (CH) 

In many cases the coefficients Cn are not of the same size, but instead they present a noticeable 
growth with the increasing perturbative order. Since the CH model relies on a specific form 


for the perturbative expansion, Eq. (4.0.1), results are not invariant under a rescaling of the 


expansion parameter from as to cts/A. Ref. [23] presents a slightly modihed version of the CH 
model. In this modified model, denoted as CH, we first rewrite the perturbative expansion of 


eq. (4.0.1) in the form 




n=l 




n=l 


with 


(n — 1)! 

(n-1)! 


(4.2.6) 


(4.2.7) 


Then we use the coefficients bn as the Cn coefficients in the CH model. The probability density 
becomes therefore 


f{Ak\bi,...,bk) ^ 


ric 


nc + lj 2kl{as/X)^+%) 
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fc!(aPA)'“+T(fe) 

lAfcl 


ric+l 


if lAfcl <fc!(^)"+'6, 
if lAfcl 6, 


(k) 


and the credibility interval 

\ k-\-l 


df = 


k\ (f)' 


6(/.) 


if 


P% < 


^ k!(f)^+H(k)[(n, + l)(l-p%)](-^/^^^ if p%> 


ric+l 


ric+l 


(k) 

(4.2.8) 


(4.2.9) 


There are several ways to determine the optimal value for the rescaling factor A. Ref. 
determines A empirically by observing how the model succeeds in predicting theoretical uncer¬ 
tainties for observables for which higher order perturbative computations are known. Another 


way of determining A is proposed in 24 , where the best A is thought to be the one that makes 


all the expansion coefficients closest to the same size. 


4.3 David-Passarino model 

Reference predicts higher orders using the concept of series acceleration. In this case they 
use of a collection of series transformations in order to improve the convergence of a series. Even 
if the original series may be (and in the QED/QCD case is) divergent, the helpful property of 
sequence transformations is that they provide a result that can be interpreted as the analytic 
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extension of the sum of the series. Therefore, the difference between the all-order sum and the 
partial sum up to the last known order is what we call the theoretical uncertainty. 

The hypotheses on which this model is based are related to the analyticity of the series. 
In fact, there are an inhnite number of functions that have the same asymptotic expansion, 
therefore we need to assume that 

1. The analyticity domain is sufficiently large 

2. After a certain perturbative order, there is an upper bound on the remainder of the series. 

The general applicability of sequence transformations to QFT theories is widely discussed in 
Ref. [^. We briefly review the most important sequence transformations that are used in the 
model. 

4.3.1 The Levin r-transform 

The first kind of sequence transformation that one can consider is the so-called Levin r- 
transform. If the partial sum has the form 


Sn = '^7i 
i=0 


we can define the r-transform as 

"vfc 


rm = 


Ei=io (ri, k, i, 13) Sn+i 
Ei=io L /3) 


Tk = T^ = r^( 0 ) 


(4.3.1) 


(4.3.2) 


where io = n — 1 and 




V ^ 1 


(4.3.3) 


where (z)a = r( 2 ; -|- a)/r( 2 ;) is the Pochhammer symbol. A represents the forward-difference 
operator: A5„ = 

The key to estimate the first unknown coefficient is to Taylor expand r^. Suppose that 
Si, ... , Sk are known. One then computes 


Tk-Sk = 7fc+i + O 


(4.3.4) 


and 7 fc+i is the prediction for 7 ^+ 1 . If the number of known values k is very small, then the 
prediction is not expected to be reliable. However, if we apply 

T 2 - S 2 = {^ 2 / 11 ) z^ + o {z^) (4.3.5) 


to the Higgs series, we end up with a result that has the correct sign and order of magnitude. 

The procedure that allows us to improve the convergence of the series is the following: 
suppose that we want to apply the Levin r-transform t°(/ 3) to the Higgs series. 


1. First use the first 3 terms in 


(|4.3.1 ), with 73 = 7 f(/r = mn), and derive 74 : 


74 


= 3 


7172 L 




/i Id 


I2 + 7/3 + /32 


(4.3.6) 
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2 . Compute S 4 assuming 74 = 74 . 

3. Derive 


75 = 


^74 
7i 72 73 


(120 + 72/3 + 15/3^ + ^^) \ 


(4.3.7) 


where 


^ = 47^73 (6 + II/3 + 6/32 + /33) -6717I (24 + 26^ + 9/32 + /33) 

+ 47i72 74 (60 + 470^ + 12/32 + /33). 


(4.3.8) 


4. Calculate ^5 assuming 75 = 75 . 

5. Repeat steps 1 to 4 until ra, ..., rg are computed. 

6 . Compare the S 3 , . ■ ■ , Sq with the T 3 , . .., tq. 

7. Repeat again steps 1-6 for 73 = 73 + A 73 and 73 = 73 — A 73 , always taken at // = uih- 

4.3.2 The Weniger 5-transform 

A second transform that is considered in is the 6 -transform introduced by Weniger: 


where 


i ('k\ if3 + 1 


W^ik,i,P) = {-l) 


+ / (/3 + A:)fc_i 7i+i++i’ 

Of course, we can follow the same exact recipe described above, and we get 


and 


75 = 


74 = 


74 


73 


37172 


(47173-72) 


10 7172 73 


(72 73 - 9 71 73 + 18 71 72 74) • 


(4.3.9) 

(4.3.10) 

(4.3.11) 

(4.3.12) 










Chapter 5 

Phenomenological applications 


In this Chapter we turn our attention to possible phenomenological applications of what we 
have learned in the previous Chapters. In particular, our first concern will be to determine the 
perturbative order at which the partial sum start to diverge. We know that at some point the 
perturbative expansion of any observable in QCD will deviate from its real value, but we do not 
know when. If, for example, a series started to diverge at the fourth perturbative order, then 
we should start to worry that our theoretical predictions might be meaningless. 

After studying the divergent behaviour of a series, one can give an estimate to its asymptotic 
value, via the Borel method illustrated in Section 1.2.3 Hence, the theoretical uncertainty on 


the perturbative expansion is defined as the difference between the asymptotic value and the 
truncated series. 

Since we have seen three sources of divergence but the instanton impact is by all means 
negligible, we will first study the effect of the Landau pole divergence and renormalon divergence 
assuming that only one source at a time is dominant. Then we will combine both sources to 
obtain the most complete estimate for the theoretical uncertainty. The combination is important 
because we know that the soft approximation is dominant at low orders, but we also know that 
the renormalon divergence will occur sooner. 

We will apply our models to two of the most relevant processes at the LHC: Higgs production 
and ti production. The former is known exactly up to N^LO 27 and its soft approximation is 
known to work very well [l^. The highest PDF order available is NNLO, so one might ask if it 


is sensible to convolute the N^LO partonic cross section with the NNLO parton distributions. 
The answer is yes, because the impact of theoretical uncertainties on PDFs is negligible for 
the Higgs cross section. For ti production, the perturbative series is known up to NNLO and 
the soft approximation has been studied in 


28 


Finally, once we have determined our theoretical uncertainty, we will have to check its con¬ 
sistency. We do so first by pretending not to know the N^LO and trying to predict it. Then, we 
compare our method with the other available methods mentioned in Chapter 


5.1 Landau pole divergence only 

The first case we consider is a model where only the Landau pole divergence is present. In this 
scenario, the all-order behaviour of the cross section is dictated by the soft approximation. In 
other words, if we take the resummed cross section in Mellin space, expand it in powers of as 
and then transform back in physical space term by term, we obtain the all-order cross section. 


This cross section has a divergent perturbative expansion like we saw in Section 2.2 
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Since the resummed cross section has the form of an exponentiation, it is easier to work 
in terms of the physical anomalous dimension, defined in Eq. (2.2.2). The physical anomalous 


dimension has the same divergent behaviour than the cross section. In particular, the order at 
which the series starts to diverge is independent of the logarithmic accuracy, so that we can 
study the behaviour of 7 ll; which is much simpler than the higher logarithmic orders. We recall 


that the series we want to study is Eq. (2.3.4) 


OO 

R{as{Q^)^x) = V AW(1) [-/3oas(Q2(l - x))]-+\ 


n=0 


(5.1.1) 


where Q can be taken as the hard scale of the process (the Higgs or twice the top mass), and x 
is what we usually call r: 

(5.1.2) 


r = 


the ratio between the invariant mass of the final state and the hadronic center-of-mass energy. 

Once that we have determined the trend of the perturbative series, our theoretical uncertainty 
is simply the difference between the resummed result and the truncated series. This time it is 
more advisable to use the most accurate result (N^LL), since while the divergent behaviour 
is independent of the logarithmic accuracy, the value of the resummed cross section strongly 
depends on it, so we cannot make any semplifications. 

Finally, |12| and 28 find that a better approximation of the fixed order result is obtained 


if we introduce certain subdominant logarithmic contributions. For the sake of clarity, we will 
not include those contributions in the soft-approximate series. We will use what is called Wsoft 
in 12 , that is to say the simplest approximation based on )D^°®(Ai) of Eq. (2.4.9). 


5.2 Renormalons only 

Now let us suppose that the dominant source of divergence is the renormalon divergence. There 
are two models that we can build based on renormalons. In the first one, which we shall call 
Naive Renormalon Model (NRM), the cross section exhibits the typical factorial divergence as 
it is. 

The second model, the Exponent Renormalon Model (ERM), is built starting from the as¬ 
sumption that the renormalon divergence can occur inside the Sudakov exponent. The physical 
interpretation behind this model is that we are considering a class of corrections not to the 
whole process, but only to the soft gluons that are emitted (e.g. fermion bubble chain along an 
emitted gluon line). 

So far, our discussion of renormalons has been made at the parton level: the divergent series 
we consider is the partonic cross section in N space. In Section |5.2.3| we will show how we can 
extend the model to hadron level results (much easier to work with) without any trouble. The 
tacit assumption is that the PDFs, which are a perturbative expansion as well, don’t present 
any divergent feature, or their divergence is subdominant with respect to that of the partonic 
cross section. 

From now on, we will suppose that the first perturbative orders of an hadronic inclusive cross 
section are known, so that the truncated N^LO cross section is given by 


'^N'=LO ~ *^0 
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n+1 
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(5.2.1) 











5.2. Renormalons only 


37 


5.2.1 Naive Renormalon Model 


Recall that when renormalons are dominant, then the Borel transform of the perturbative series 
has the form 

°° K 

B[fl](t<)= ^ ( 5 . 2 . 2 ) 


m=—oo 


u ’ 
m 


where u = /3ot- Renormalons are poles located at u = m for every integer m, but the values 
of the residues Km are unknown. There are some ways to determine these residues from first 
principles, using the bubble diagram chains, but those methods rely on strong approximations. 
The exact determination of the residues would require an all-order calculation, which is precisely 
what we are trying to avoid. Instead, our approach consists in exploiting what we know about 
the series (the hrst known coefficients), to extrapolate the first and most important residues by 
comparison. 

If we know the cross section Eq. (5.2.1) up to N^LO, and therefore we possess three pieces 
of information beyond LO, we can assume that the Borel transform of the series only has three 
poles: 


B[R]{u) = 
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m 3 


(5.2.3) 


1 — u/mi l — u/m 2 1 —u/ma’ 

where irin are the integers representing the position of the renormalon poles in the u axis. If 
this is the Borel transform, then the series has the form 
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with the coefficients given by 
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(5.2.5) 


Now we ask that the NRM predict the first known orders exacly and solve the linear system for 
77 G [0, 2] to determine the residues: 


Cq — Kmi T Km2 T 


‘■m 3 
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(5.2.6) 
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where Cj are the coefficients of Eq. (5.2.1). Once that we have the residues, we know the 
perturbative series at all orders in and we can easily see where the series begins to diverge. 

In order to evaluate its asymptotic value, recall that in Section 1.2.3 we defined the Borel 
integral as (3.1.12), so that in our case we can assign the following value to the sum of the series: 

/ 3 
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5.2.2 Exponent Renormalon Model 

The ERM works exactly as the NRM, with the only difference that the renormalon series expo¬ 
nentiated. Let us take as an example the Higgs cross section. Recall that the resummed cross 
section has the form: 

= aogoias)expE{as). (5.2.8) 

In the ERM, the exponent has the typical renormalon form 

OO 

E(a,) = (5.2.9) 

n=0 


with rn given by Eq. (5.2.5). To extract the values of the residues, we expand Eq.( 5.2.8) in 
powers of as'- 
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(5.2.10) 


and then we confront order by order with the exact cross section to extrapolate the residues 


1 and Rims • 


Once again, using the residues we are able to know the series at all orders. The asymptotic 
sum of the series is nothing but the exponential of the Borel integral: 
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(5.2.11) 


5.2.3 Partonic vs hadronic cross section 

We have formulated our models at the hadronic level, but to be more precise we should have 
applied our models to the partonic cross section in N (or z) space. In fact, for every N there 
should be a renormalon series with residues Ki{N). We will now show that we can apply both 
renormalon models directly to the hadronic cross section, without going through N space. The 
basic idea is that once we perform the inverse Mellin integral to obtain the hadronic cross section. 


the ^dependence is converted to a r-dependence. But r (5.1.2) is fixed by the invariant mass 


of the final state and by the center-of-mass energy of the collision. Hence, the residues are pure 
numbers. 

Let us consider the NRM with only one renormalon for simplicity. The generalization to the 
case of multiple renormalons and to the ERM is straightforward. Remember that the inclusive 
hadronic cross section, differential only in M, can be written as the inverse Mellin transform of 
the cross section in N space: 
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where the hadronic cross section in N space cr{N) is given by 

'f3o 


K{N)Y, 


n=0 
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a{N) = ao{as)^{N) 

Now, once we perform the Mellin inversion integral we obtain 

OO 

a{T,M'^) = fTo(Q;s) ^ 
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Therefore, the residue at hadronic level can be defined as 
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(5.2.14) 


(5.2.15) 


The same exact reasoning holds for the physical space (although the coefficient functions are 
distributions) and, with small modifications, for the ERM. 


5.3 Renormalons and Landau pole divergence 

Since the soft approximation is particularly good at all the known orders, we can make the 
assumption that the series is dominated by soft terms at very low perturbative orders, and 
the renormalon effects start to kick in later. In this scenario, the renormalons would cover the 
discrepancy between the soft approximation and the exact result at low orders. Once again, we 
divide the case in which the renormalons are outside or inside the Sudakov exponent. 


5.3.1 Naive Soft-Renormalon Model 


Let us consider a situation where the hadronic cross section is simply the sum of a soft and a 
renormalon contribution: 
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(5.3.2) 


is the soft approximation to the exact cross section based on of Eq. (2.4.9). More 

specifically, considering the N^LO cross section of Eq. (5.2.1), we can separate the soft and the 
renormalon part: 
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n=0 


We now compare the renormalon coefficients of Eq. (5.2.5) with the renormalon coefficients 
of Eq. (5.3.3), that we can determine recursively from a and using 
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For what concerns the asymptotic value, the sum of the soft part of the series is the resummed 


NNLL result computed with the Borel prescription described in Section 2.3, Note that since 
we are interested in the soft contribution only we shall not use the NNLL+NNLO result like in 


Eq. (2.1.19). We denote as resummed Wsoft cross section (2.4.9) without fixed-order 

matching. The sum of the renormalon part, instead, is the usual Borel integral. Therefore, the 
sum of the whole series in the Naive Soft Renormalon Model is 
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5.3.2 Exponent Soft-Renormalon Model 

The equivalent of the ERM that fully includes soft contributions is a model where the cross 
section has the form 
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where we used the renormalon exponent of Eq. (5.2.9) and < 7 ®°^ of Eq. (5.3.2). To compute 


the residues, we need to compare the fixed order results with the perturbative expansion of 


Eq. (5.3.6), which is 
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Then, the asymptotic sum of the series is 
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5.4 Divergence point n 

In this Section we analyze the trend of the partial sum in each of the models, in order to see at 
which perturbative order the divergence occurs. We will see that the Landau pole divergence 
occurs much later than the renormalon divergence. We will also show that the divergence point 
in any of the renormalon models only depends on the location of the leading renormalon pole. 
It goes without saying that in the models with both sources of divergence, the renormalon 
divergence prevales. 
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5.4.1 Landau pole divergence only 


Let us recall the form of 7 ll Eq. (5.1.1) and write its partial sum Rj^: 


Rk = '£^n = Yl [-PoasiQ^l - X))] 


n+1 


n=0 


n=0 


where A{z) = l/r(2:). This series represents the physical anomalous dimension at LL. The 
physical anomalous dimension is related to the exponent of the resummed cross section, and it 
is obvious that if a series has the form an oc ex.p{En), then an starts to diverge when En starts 
to diverge. 


In Fig. 5.1 the trend of the partial sum Rk and of the sequence of contributions Sn is shown. 
We see that the series converges rapidly to a plateau (the quantity that is summed at each 
perturbative order is very small), and then the partial contributions start to gain a noticeable 
value at very high perturbative orders {k ~ 70). One usually defines the perturbative order n 
at which the series starts to diverge by imposing that 


ISn—1| 


for n > n. 


(5.4.1) 


With this definition, we find that h = 33. Note that this does not imply that the series is already 
divergent at h. It only means that from that point on the contributions that we are adding start 
to grow. However, the macroscopic effect is only visible at n ~ 70. 


The result only depends on x, which we take as r of Eq. (5.1.2). The x value corresponding 
to the worst behaviour of 7 ll is x ~ 1. However, for the Higgs production at 13 TeV center- 
of-mass energy x ~ 10“^. Therefore, we conclude that there is no reason to worry about the 
Landau pole divergence since, for processes far away from threshold, the divergence kicks in at 
very high perturbative orders. 


5.4.2 Renormalon divergence only 

For both models that involve the renormalon divergence, the order at which the series starts to 
diverge h depends on the location of the leading renormalon. Before plotting the outcome of 
the models, we will justify this statement analytically. In our models multiple poles are present 
but, as we are going to show, the dominant pole is the one that is closer to 0. Therefore, let us 
suppose that the renormalon series has only one pole in its Borel transform and hence only one 
residue. The partial sum in this scenario is 

n 

Rn = KY,{aPofk\a^+\ (5.4.2) 

A:=0 

where a = I/m is the inverse of the location of the pole in the u axis. To prove our point, we 
shall compute the increment and study its trend with respect to n: 

Rn+i - Rn = Kaf3oal {a/3oas)'^T{n + 2). (5.4.3) 

If we continue this function analytically for real n and find a minimum at no, it means that 
the contributions that we are adding to the partial sum are decreasing before no, and start 
increasing afterwards. To look for a minimum, we compute the derivative of the increment with 
respect to n and put it to 0. 


^ [(a/3oas)"' F(n + 2)] = 0 log(a,doas) + V’(^ + 2) = 0, 


(5.4.4) 
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Figure 5.1: Trend of the partial sum (top, blue) and of the sequence Sn (bottom, red), 

TnP" 

supposing that only the Landau pole divergence exists. Results are obtained using x = 
with s = 13 TeV, and are almost identical if we use the ti invariant mass instead of the Higgs 


mass. 
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where ipi^) is the digamma function, the logarithmic derivative of r(z). We could solve this equa¬ 
tion numerically, but a very interesting formula comes out if we use the Stirling approximation, 
supposing that n is big, as can be checked a posteriori. In fact, if we define t = n + 1, 


'ip{t + i) = iogr(t +1) ~ ^log 


Now, since we are considering large t 


(-Xv^t 

d 

\e) 

dt 


tlogt-t+ ^logt 






and therefore we have 


1 


-g 1) ~ log { 1 -g — 


The equation which determines the minimum becomes therefore 


log ( 1 + ^ j = log ^ 


aPoa, 


1 


n = 


aPoOs 


(5.4.5) 


(5.4.6) 


(5.4.7) 


(5.4.8) 


where m denotes the location of the pole in the Borel transform. Numerical solutions of 


Eq. (5.4.4) for integer values of 1/a are in good agreement with this approximation. 

We have thus proven that the divergence point is directly proportional to the position of the 
leading renormalon. 


There are arguments 
like Drell-Yan, Higgs and tt production 


14 in favour of the hypothesis that, for processes 
the leading renormalon is located at m = 1. These 


arguments are based on the fact that a previously supposed pole at m = 1/2 is in reality 
canceled at all orders. Hence, we conclude that the renormalon divergence occurs much before 
than the Landau pole divergence, more precisely at 


n ~ 13. 


(5.4.9) 


Finally, we note that Eq. (5.4.4) has the same solution if we consider the ERM instead of the 


NRM, as one could easily expect. 

One last observation can be made about the instanton divergence, justifying what has been 
said in Section 3.2.1 Since the location of the leading instanton pole in the Borel plane is at 
t = 47r, the instanton divergence occurs at h ~ 110, therefore much later than the Landau pole 
divergence. 


5.5 Theoretical uncertainties 


This Section contains our predictions for the theoretical uncertainties on the perturbative series 
in QCD, based on the models that we introduced above. In each of the models, the difference 
between the all-order result and the last known order can be used as an estimate for the theo¬ 
retical uncertainty. However, there is a certain arbitrariness in determining the asymptotic sum 
of the series, i.e. the freedom in choosing the renormalon poles whose residues are extrapolated 
by comparison from the exact known orders. 

We will consider in particular the process of Higgs production in gluon fusion and, in Sec¬ 
tion 5.7, ti production in gluon fusion. Throughout this Section, for phenomenology we will 














5.5. Theoretical uncertainties 


44 


Table 5.1: Values of the Borel integral in the NRM for the Higgs cross section at 13 TeV, pdf set 
MSTW2008nnlo68cl. Results up to NNLO obtained with ihixs, the N^LO is taken from 


27 


Pole choice (mi m2 m3) 

Borel integral [pb] 

(1 2 3) 

44.5645 

(1 3 4) 

44.419 

(2 3 4) 

43.8916 

(-2 -3 -4) 

42.4534 

(-1 -3 -4) 

42.1961 

(-1 -2 - 3 ) 

42.0315 


use the code ihixs 1.4 to obtain the inclusive Higgs cross section up to NNLO at a scale 
fJ-R = ^J‘F = TRiH- We will then add the exact N^LO result that can be found on 


27. This 


last result is obtained using the PDF set MSTW2008nnlo68cl, so we will consistently use that 
same PDF set as an input for the code. Since we are trying to neglect all the PDF perturbative 
effects, we will always use the most accurate set available, i.e. NNLO. The impact of theoretical 
uncertainties on PDFs on the Higgs cross section is negligible |24| . 

We are going to need the soft-gluon approximation to the Higgs cross section A^-soft, de- 
To compute it, we use ggHiggs 2.1 


scribed in Section 2.4 
result we use TROLL 
described in 


30 , while to compute the resummed 


31 with the V^-soft prescription. Remember that the soft approximation 


12 also contains some sub-leading contributions, but for the sake of clarity we are 


not going to take them into account. We will use the simpler approximation based on 


as in Eq. (2.4.9) 


For tt production, we will use Top++ 2.0 32 with NNPDF3.0 PDFs. The V'-soft approxi¬ 


mation, discussed in Ref. will be extracted from a private code by Claudio Muselli. 

5.5.1 Choice of the poles 

Like we said before, for the Higgs case only three orders are known beyond the LO, so we have 
a system with 3 equations and 3 unknowns residues. The location of the poles of Eq. (5.2.5), 


however, is free. For example, an obvious choice would be mi = 1, m 2 = 2 and m 3 = 3, but then 
our model would not predict the behaviour induced by the UV renormalons, which are located 
in the negative u axis. 

If we compute the asymptotic sum of the series for different choices of the poles, we notice 
that the maximum value Emax is obtained for (mi m 2 m 3 ) = (l 2 3), and the minimum 
Emin for (—1 —2 —3), all the other values being included inside that interval. This is not 

easy to prove, but one can solve the system (5.2.6) for generic m, and see that the results of the 


Borel integral Eq. (5.2.7) obey to this rule. 

In Table [RT} as an example, the asymptotic sum of the Higgs cross section is computed using 


different pole choices in the NRM. 


5.5.2 Trend of the perturbative expansion 


In Fig. 5.2 the partial sum of the series an in the NRM is plotted as a function of the perturbative 


order n. The different colors represent different choices of the poles, in particular we plot the 
one that maximizes and the one that minimizes the asymptotic sum of the series. We note that 
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Figure 5.2: Trend of the partial sum as function of the perturbative order for the Higgs cross 
section at 13 TeV center-of-mass energy in the NRM. Results up to NNLO are obtained with 
ihixs, the N^LO is taken from 
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PDF set MSTW2008nnlo68cl. 


the series exhibits a fixed-sign divergence when IR renormalons are present, while it exhibits an 
alternating-sign divergence when UV renormalons are present, like we studied in Chapter 

We have chosen the NRM as en example: plots like Fig. 5.2 are similar in all the models. 
What really changes between different models is the estimate for the theoretical uncertainty on 
the last known order. In the next Section, we will compare the theoretical uncertainty estimates 
between different models. 


5.5.3 Comparison between our models 

Let us summarize how we can study the higher-order behaviour of a perturbative expansion in 
QCD. In the NRM, we assume that the series has a simple factorial divergence. In the ERM, the 
factorial behaviour is exponentiated. The NSRM and ESRM are the extensions of respectively 
the NRM and the ERM when we take into account the soft-gluon approximation. 

The trend of the partial sum does not depend strongly on the model: the leading renormalon 
pole is always located at |ti| = 1 and therefore the divergence occurs always at h ~ 13, like we 


saw in Section 5.4.2 In fact, even when the Landau pole divergence is present together with the 


renormalon divergence, the latter is dominant. 

The freedom to choose the poles is also common to each of our models. In particular, 
it is always true that the maximum and minimum values for the asymptotic sum of the series 
correspond to the choice of the poles (1, 2,3) and (—1, —2, —3), although maximum and minimum 
can be exchanged. What matters is that any other pole choice leads to an asymptotic sum 
included in that range. 

Once that we have the truncated sum (the N^LO) and the asymptotic value S, we say 
that the difference between the two is our theoretical uncertainty. However, we do not have a 
single asymptotic value, but a range of values [Emin, Smax]- Our uncertainty interval A, then, is 
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Figure 5.3: Higgs cross section: intervals for [Smax, Flmin], the maximum and minimum value for 
the asymptotic sum of the series, together with the exact N^LO. Results from different models 
are compared. We used ihixs and [27] , finite top mass and finite bottom mass effects are 
included and the pdf set is MSTW2008nnlo68cl. 
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^ [tTfc , Sjnax. 
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if (Tfc > Sjuax 
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(5.5.1) 


Now, what really makes our models different from one another is precisely the value of the 
predicted S m ax and Smin- In Fig. |5.3[ these asymptotic values are computed in each of our 


models and compared. At a first glance, we notice that the Exponent Models tend to give a 
wider spread of possible values for the asymptotic sum S, and therefore will have a broader 
theoretical uncertainty band than the Naive Models. 


5.6 Accuracy of the models 

We have introduced 4 different models for the theoretical uncertainty on the Higgs cross section, 
but their results don’t seem to be in good agreement. In order to determine which of them give 
accurate predictions we perform two test. First we see which models describe the uncertainty 
on the NNLO correctly, then we compare our estimates with the estimates predicted by all the 
other models described in Chapter 

5.6.1 Accuracy in predicting the N^LO 

Let us pretend that we do not to know the N^LO nor the resummed N^LL cross section. We can 
apply our methods to the truncated partial sum up to NNLO and see if the uncertainty bands 
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Figure 5.4: Theoretical uncertainty bands on the NNLO Higgs cross section (blue lines), com¬ 
puted pretending not to know the N^LO. The exact N^LO is displayed as a dashed green line. 


include the true N^LO or not. 

If we know the cross section up to NNLO, there are only two residues that we can fit, and 
therefore all the models suffer from lack of precision. The Smax and Emin values are obtained us¬ 


ing (1, 2) and (—1, —2) as poles. In Figure 5.4, results are plotted for the theoretical uncertainties 


on the NNLO. The plot also displays the actual value of the N^LO. 

We notice that the only model whose uncertainty band contains the N^LO cross section is 
the Exponent Renormalon Model. However, the uncertainty band in the ERM is very large. 
The NRM and NSRM are not far from the true value, while the ESRM seems to completely 
miss it. 


Eurthermore, if we compare Eig. |5.4| and Eig. 5.3, we notice that the ERM and the ESRM 


have completely opposite predictions. In fact, if we include the N^LO the ERM predicts an 
asymptotic sum which is smaller than the N^LO itself. Without N^LO, instead, the uncertainty 
band is on the upper side of the NNLO. The opposite thing happens for the ESRM. We have 
checked that, due to the exponential, the outcomes of the Exponent Models are highly unstable 
with minimal variations of the initial values of the cross section. Eor example, if we use the cross 
section up to NNLO computed by ggHiggs, we obtain a very short uncertainty band and the 
approximate N^LO lies perfectly inside that bancQ 

We conclude that the Naive Models represent our best models for THU. They are both quite 
accurate at NNLO and stable with small variations of the initial conditions. Eurthermore, if we 


compare Eig. 5.3 with Eig. 5.4, we note that the uncertainty band on the N^LO is smaller than 
the one on the NNLO, as one would expect from a convergent method. In particular, the NSRM 
predicts a positive uncertainty band, which we expect from the monotonicity of the perturbative 
cross section. 


^The NNLO computed with ggHiggs differs from the NNLO computed with ihixs by some finite top mass 


effects. 
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Figure 5.5: Various models for the theoretical uncertainty on the N^LO are considered. In 
red the Scale Variation uncertainty is displayed, in violet the Cacciari-Houdeau modified model 
CH with 68% C.I., in green the David-Passarino uncertainty and in blue our models: the Naive 
Renormalon Model and the Naive Soft Renormalon Model. 


5.6.2 Consistency with other models 


After checking the accuracy of our models at NNLO, we have ruled out the Exponent models 
since their predictions don’t seem to be reliable. We would like now to compare the estimates for 
the THU in the NRM and in the NSRM with the other known models described in Chapter 
There, we saw that the conventional way to give a theoretical error band is to vary the 
renormalization and factorization scale fiR and fip around a central value (usually = fip = 
mu)- This model is called Scale Variation. The Cacciari-Houdeau model [^, which we shall 
use in its modified version for hadronic observable [23] , is instead based on a Bayesian approach 
and gives a 68% confidence level interval for the theoretical uncertainty on higher orders. The 
David-Passarino model 


of the serie^ 


25 makes use of sequence transformations to improve the convergence 


Fig. 5.5 shows the comparison between all the models. The NRM appears to be consistent 
with the Scale Variation, while the NSRM has a positive uncertainty band like the DP model, 
but its width is much shorter. 


^The uncertainty given in |25j is valid for the Higgs cross section at 8 TeV, here we have simply rescaled the 
uncertainty band for the 13 TeV cross section. 
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5.7 tt production 


What has been said so far is valid for the cross section of Higgs production in gluon fusion at 
13 TeV. We turn now our attention to another relevant process at the LHC: tt production in 
gluon fusion. There are some important differences between the tt and the Higgs case. First, 
the maximum available order for ti is NNLO, which means that we can at best extrapolate 
two residues. Then, the THU on the PDFs has a larger impact on tt than on the Higgs cross 
section 24 . Finally, the soft approximation for tt is slightly worse than the Higgs case in 
predicting the first known orders |28| . 

For what concerns our models, we can apply them to the tt hadronic cross section with¬ 
out significant modifications, apart for the ERM. In fact, the resummed cross section has the 
following form 




5i(as)exp[G'i(iV)]-hO 
1 = 1.8 ^ 


(5.7.1) 


where we can see the separation between the singlet and the octet components. The constants 
5i(cts) can be found in (^ . Accordingly, we need to modify Equation (5.2.9) by introducing 
two renormalon exponents: 


Ei{as) = 


OO 

n=0 


a 




Esias) = 


OO 

n=0 


a 


n+l 


(5.7.2) 


Since we have the sum of two exponents, we extrapolate one residues from the first exponent 
and one from the second exponent. However, we don’t have to worry about this in the ESRM 
since the model can be formulated with Eq. (5.3.3), where already includes the sum of 


the singlet and the octet contributions. Therefore, we are left with only one exponent for the 
renormalon part of the series. 

We have tested the accuracy of our models at NLO, the penultimate available order. Results 


are shown in Eig. 5.6 (the equivalent of Fig. 5.4). Once again, we decide to trust the Naive 
Models. In particular, the NSRM, our preferred choice for the Higgs cross section, is the most 
accurate in predicting the NNLO. 


Finally, Fig. 5.7 shows the theoretical uncertainty bands on the NNLO. Comparing Fig. |5.7| 
with Fig. |5.6[ we notice that once again the ERM model is unstable and changes drastically 
when a new order is included. The Naive Models have a behaviour which is very similar to the 
Higgs case. However, the NRM uncertainty band unexpectedly grows when a perturbative order 
is added. This does not necessarily mean that the NRM model does not converge (i.e. predicts 
higher and higher uncertainties). We recall that at NLO we only have one residue to determine, 
so we can expect an underestimation of the uncertainty since such little information is available. 

We conclude that even for the tt case our preferred model is the most complete, accurate 
and convergent one: the Naive Soft Renormalon Model. 
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Figure 5.6: Theoretical uncertainty bands on the NLO ti cross section (blue lines), computed pre¬ 
tending not to know the NNLO. The exact NNLO is displayed as a dashed green line. The fixed- 
order results and the resummed result are provided by Top++, the PDF set used is NNPDF3.0 
with as = 0.118. 



Figure 5.7; Theoretical uncertainties on NNLO ti cross section (blue lines), predicted with each 
of our models. The NNLO is the last available known order for tt production. The settings are 
the same as the Figure above. 










































Conclusions 


In this thesis, we have studied the divergent behaviour of the perturbative series in QCD, taking 
as an example the Higgs cross section and the ti cross section. We have seen that there are three 
sources of divergence: the Landau pole, the renormalons and the instantons. After introducing 
each of the sources, we have studied their phenomenological impact on the all-order cross section. 

In some theories, like the double-well potential in QM, the instanton divergence can com¬ 
pletely cancel the divergence of the perturbative expansion. However, this is not the case of 
QCD because of the IR renormalons. We found that the impact of the instanton divergence is 
totally negligible, as the partial sum an would start to increase signihcantly at h ~ 110. 

The Landau pole divergence, related to the soft-gluon resummation, occurs at high orders 
too (n ~ 33). However, even if we don’t have to worry about the divergence, we can construct 
models that combine soft contributions with renormalon contributions. We do so because we 
know that, accidentally, the expansion of the resummed cross section in the soft limit provides 
a good approximation of the exact known orders. 

Finally, the dominant source of divergence in QCD is the renormalon divergence. Renor¬ 
malons are poles in the Borel plane, which introduce an ambiguity in the definition of the 
asymptotic sum of the series. The divergence point for renormalons depends on the location of 
the leading pole, and in our case it turns out to be h ~ 13. 

After studying the divergent behaviour of the partial sum, we computed the asymptotic 
value of the sum of the series S via the Borel integral method. We defined our estimate for 
the theoretical uncertainty as the difference between the last known order and S. However, the 
renormalon models allow a certain freedom in predicting S. In fact, the residues at the poles of 
the Borel transform cannot be computed theoretically and need to be extracted by comparison 
with the exact fixed-order results. The value of S depends on what residues we decide to extract. 
The Borel transform has poles for each integer value of the real axis, but we can determine only 
3 (Higgs) or 2 (tt) residues. Therefore, combining the range of possible values for S with their 
distance to the last known order, we can give an estimate to the theoretical uncertainty. 

In particular, we have formulated 4 different models corresponding to 4 different ways of 
seeing the renormalon effect. First, we can simply assume that the perturbative expansion of 
the Higgs cross section presents the renormalon factorial divergence as it is (Naive Renormalon 
Model). Otherwise, we can assume that the renormalon corrections (basically identified with 
chains of bubble diagrams) are applied on the emitted soft gluons, and therefore the factorial 
behaviour is exponentiated (Exponent Renormalon Model). Finally, for both models we can 
consider the extension that takes into account the soft-gluon approximation. In other words, 
since we know that the soft-gluon approximation dominates at low orders, we ask that the 
renormalon divergence should cover the difference between the exact fixed-order series and its 
soft-gluon approximation. Those models are called Naive Soft Renormalon Model and Exponent 
Soft Renormalon Model. 
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We have compared the theoretical uncertainty predictions from our 4 models and we have 
found out that the uncertainty bands in the Exponent Models are bigger than in the Naive 
Models. Then we have investigated the accuracy of our models by computing the uncertainty 
bands on the penultimate known order and checking if they predict the last known order correctly. 
In this investigation, we have noticed that the Exponent models are unstable for small variation 
of the initial parameters. So we have ruled them out and we have compared the NRM and 
the NSRM with the other already known models for theoretical uncertainties. Our models are 
compatible in magnitude with the conventional Scale Variation model. Actually, the NSRM, 
unlike the NRM and the Scale Variation, predicts a positive-only uncertainty, which one could 
expect since the perturbative series for the Higgs and the ti cross section appear to be monotonic. 

Finally, we have compared theoretical uncertainties on the two different processes, finding 
similar results. Our preferred choice is the Naive Soft Renormalon Model, since 

• It is the most accurate in predicting the last known order 

• The uncertainty gets smaller when perturbative orders are added 

• It is stable for small variations of the initial conditions 

• It predicts a positive-only uncertainty, compatible with a monotonic partial sum. 

In the NSRM, the Higgs and ti cross sections at the LHC at 13 TeV become 


app^H = 43.14+g^^"^^" pb. 

(5.7.3) 

app^u = 720.195+2^^5% p^. 

(5.7.4) 
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